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We study the magnetic properties of spherical Co clusters with diameters between 0.8 nm and 
5.4 nm (25 to 7500 atoms) prepared by sequential sputtering of Co and AI2O3. The particle size 
distribution has been determined from the equilibrium susceptibility and magnetization data and 
it is compared to previous structural characterizations. The distribution of activation energies was 
independently obtained from a scaling plot of the ac susceptibility. Combining these two distribu- 
tions we have accurately determined the effective anisotropy constant Keff- We find that K^ff is 
enhanced with respect to the bulk value and that it is dominated by a strong anisotropy induced at 
the surface of the clusters. Interactions between the magnetic moments of adjacent layers are shown 
to increase the effective activation energy barrier for the reversal of the magnetic moments. Finally, 
this reversal is shown to proceed classically down to the lowest temperature investigated (1.8 K). 
PACS:75.50.Tt,75.70.-i,75.40.Gb,75.30.Pd 
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I. INTRODUCTION 



Single domain magnetic particles are attractive for applications in data storage. Their properties differ from those 
of the bulk magnetsEl because, as the size of the particles decreases, an increasing fraction of the total magnetic atoms 
lies at the surface. The electronic raiad magnetic structure of these atoms can be modified by the smaller number of 
neighbors as compared to the bulka^cl and/or by the interaction with the surrounding atoms of the matrix where the 
particles are dispersed. For example, it was shown by Van Leeuwen et al. that the bonding of, CO at the surface of 
Ni clusters induces quenching of the magnetic moments of those atoms located at the surface.H In some other cases, 
the surface layer is oxidized and shows antiferro or spin-glass like arrangement of the magnetic moments, which also 
leads to a smaller net magnetic moment of_the particlcQ By contrast, measuring "bare" particles of Fe, Co and Ni 
produced in beams, de Heer and coworkers!] found that the net magnetic moment per atom increases as the size of 
the cluster decreases, ajmroaching the limiting value for a free atom. In addition, the net anisotropy of the particle 
exceeds the bulk value. *eL This excess was recently correlated to the augmentation of the orbital magnetic moment of 
the peripheral atoms J13'tll 

Magnetic nanoparticles are also good candidates for , t he study of quantum effects in intermediate scales between 
the microscopic and the macroscopic classical world.EHIEj In real systems however, we usually deal with macroscopic 
ensembles of particles with different sizes and shapes. The average magnetic properties of these systems come from 
intra-particle as well as interparticle phenomena, which are usually difficult to disentangle. Therefore, in this field of 
research it is desirable to obtain systems in which each of the parameters, such as the average particle size, the particle 
size distribution, the crystalline structure, and the spatial arrangement of the particles, can be varied independently 
of each other. 

We believe that the work reported here is a step forward in this direction. We present the magnetic characterization 
of a new type of systems of Co Ji^noparticles, embedded in an amorphous matrix of AI2O3, prepared by sequential 
deposition of both materials. EjIIj By varying the deposition time, the diameter of the aggregates can be controlled 
between below 1 nm and 7 nm. An important advantage of this preparation method is that it gives a rather ho- 
mogeneous dispersion of the particles in the matrix. It was also found that, for a range of thicknesses, a relatively 
ordered disposition of the partieles is obtained, in which they are arranged in layers separated from the adjacent 
ones by a controllable distance.ll3 The paper is organized as follows. In the first two sections we briefly describe the 
method employed to prepare the samples and their physical characterization. Then, we present our experimental 
results. Using the data obtained from ac magnetic susceptibility, zero-field cooled (ZFC) and field-cooled (FC) mag- 
netization measurements, and isotherms of magnetization as a function of the field, we have determined the particle 
size distribution in samples which have been prepared with different Co deposition times. We compare these results 
with available data from a previous structural characterization. This important information is then used to determine 
accurately the effective anisotropy constant and its variation with the size of the particles. We have also been able to 
separate surface anisotropy effects from the effect of the dipole-dipole interaction between the magnetic moments of 
the particles. The last section is left for the conclusions. 



II. MORPHOLOGY AND STRUCTURE OF THE SAMPLES 



Details of the sample preparation and of its structural characterization have already been reported elsewhere.l 




The Co aggregates were prepared by sputter deposition of Co atoms on a smooth alumina surface. The amount of 
deposited Co is given here by the nominal thickness tco that the deposits would have if they were homogeneous. This 
amount was measured by using energy-dispersive X-ray spectroscopy in the transmission electron microscope, and 
found to be within less than 5% of the planned dose in all cases. Clusters are formed below the percolation limit 
which appears to occur at tco = 2 nm. On top of each Co layer a new alumina layer of about 3 nm was deposited. 
Oxidized Si was used as a substrate. —A given sample is usually made by piling up a number A^ (1 to 100 for the 
samples studied here) of these layers.EEl The deposition rates of both Co and alumina were respectively 0.114A s~^ 
and O.43A s~^. It was found that the amount of Co deposited on the surface is larger than the Co mass which forms 
clusters visible by transmission electron microscopy (TEM) . The relative difference between these quantities increases 
as tco decreases. Therefore, we have in our samples non-aggregated atoms or very small clusters, which contribute to 
the magnetic signal of the samples, in addition to Co aggregates. One of the difficulties of the interpretation of the 
magnetic data is to separate these two contributions. 

Because of the upper alumina layer the aggregates show no trace of oxide even after exposure to air. The lack of oxide 
was checked by electron energy-loss spectroscopy (EELS), X-ray photoelectron spectroscopy and X-ray absorption 
spectroscopy. The atomic structure of the clustepi^was determined by extended X-ray absorption fine structure 
(EXAFS) spectroscopy and high resolution TEM.Efl It was found that the particles bond poorly to the alumina 
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matrix, and that the Co crystaUizes in the fee phase for tco < 1 nm. The presence of a fee phase in place of the. 
hep phase which is the stable phase for bulk Co is not uncommon for small particles. It was theoretically predictedtZI 
and also found experimentallycSl that the fee phase becomes more stable below some diameter which depends on the 
matrix. 

The morphology, size and spatial distribution of the aggregates were also studied using the TEM data. The 
aggregates are of approximately spherical shape (at least for tco < 1 nm) . The average diameter {D) of the particles 
increases linearly with tco- We give in Table I a list of the important parameters obtained from these experiments 
for all samples studied. Finally, the TEM pictures reveal a quasi-ordered arrangement of the Co clusterst3 that is 
induced by the topology of the layers: the clusters of a layer nucleate preferentially in the hollows left by the previous 
layer. In each layer, the average distance between the borders of adjacent clusters is of order 2 nm and approximately 
independent of tco ■ 



III. EXPERIMENTAL DETAILS 



The magnetic measurements were performed using a commercial SQUID magnetometer. The temperature range 
of the measurements was 1.8 K < T < 320 K and magnetic fields up to 5 T could be applied by means of a 
superconducting magnet. The ac susceptibility was measured by applying a small ac field (4.5 Oe) to the sample 
and using the ac detection option of the same magnetometer. The frequency lo /2t: of the ac magnetic field can be 
varied continuously between 0.01 Hz and 1.5 kHz. The samples had a rather large diamagnetic signal arising from 
the silicon substrate. This contribution was estimated independently by measuring a bare substrate and found to 
be linear in field and independent of the temperature. It was subsequently subtracted from all experimental data. 
Unless indicated otherwise, the data shown in this paper were measured on samples having more than 20 C0/AI2O3 
bi-layers in order to maximize their magnetic signals. We have checked for tco = 0-3 nm and tco = 0-7 nm that the 
variation of the magnetization and of the ac susceptibility with temperature and magnetic field is rather insensitive 
to the precise value of N, provided that N is larger than 10 layers. 



IV. RESULTS AND DISCUSSION 



A. Superparamagnetic blocking 

The magnetic dc susceptibility was measured by cooling the samples in zero field (ZFC) or in the presence of the 
measuring magnetic field (FC). Typical ZFC-FC magnetization curves are plotted in Fig. 1. At high temperatures 
the ZFC and FC curves coincide, indicating that the samples behave as superparamagnets. In this region, both curves 
follow the Curie- Weiss law C /{T — 9). The value of C increases as tco increases (see the inset of Fig. 1), as expected 
for larger clusters formation as the deposition time of Co increases. The Curie- Weiss temperature 9 is nearly zero but 
for the two samples containing the largest particles. This is brought about by the interaction between the particles, 
which we shall consider in a separate section below. At lower temperature, the two curves start to separate. The ZFC 
curve shows a maximum at a temperature Tb below which the magnetic moments are blocked in fixed directions. It is 
well known that the phenomenon of blocking is related to the magnetic anisotropy of the particles. tj The anisotropy 
favors some particular orientations of the magnetic moment, two opposite to each other in the simplest case of uniaxial 
anisotropy, which are separated by activation energy barriers U . As the temperature decreases, the number of thermal 
phonons of energy equal or larger than U xlecxeases, thus leading to an exponential increase of the time r needed to 
reverse the magnetic moment of a particleE3~0 

T = To cxp [U/kBT) (1) 

Here tq ~ 10^ — 10^^^ s is an inverse attempt- frequency, which depends on the damping of the magnetic moment 
by the phonon or the magnon baths. In this simple picture, the superparamagnetic blocking takes place when r 
equals the measurement time of each experimental point t^, thus Tb — aU/kB ln(te/'''o), where a is a constant which 
depends on the width of the particle size distribution (more details are given below). We have indeed observed that 
Tb increases with the Co deposition time, that is, with the average volume of the aggregates. Therefore, we write 
U — K^ffV, where K(,ff is an effective anisotropy constant with contributions from the intrinsic magnetocrystalline 
anisotropy of the fee Co and from other sources, such as stress indu ced anisotropy or surface induced anisotropy. The 



dependence of Kf,f / on will be considered below in section [V C , 

As expected, the blocking of the magnetic moments by the anisotropy also leads to a maximum in the temperature 
dependence of both the real and the imaginary components of the ac magnetic susceptibility. A typical experimental 
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result is shown in Fig. 2. The position of the susceptibihty peak shifts towards lower temperatures as the frequency 
of the ac magnetic field decreases since te equals 1/oj. 

Above Tb the magnetization isotherms are fully reversible because the magnetic moments are in thermal equilibrium. 
As shown in Fig. 3, the experimental data measured well above the blocking temperature of each sample collapse 
into a single curve when they are plotted as a function of H/T, indicating that the effect of the anisotropy is weak. 
Furthermore, pure Langevin curves fit the experimental data reasonably well, which shows that the size distributions 
of all these samples are narrow. Below Tb, the magnetization shows hysteresis (see Fig. 4) with both the coercive 
field He and the remanence Mr increasing as the temperature decreases (see Fig. 5). 

We plot in Fig. 6 the low-T values of Mr and of the saturation magnetization Ms as a function of tco- It is 
interesting to note that, for tco < 0.7 nm, the reduced remanence rrir = Mr/Mg is smaller than the value 1/2 
predicted by the Stoner-Wolhfarth modelnS We attribute the decrease of to a paramagnetic contribution, which 
adds to that of the blocked particles. A Curie tail shown by the saturation magnetization at the lowest temperatures 
(cf Fig. 5) is also related to this extra contribution. The excess paramagnetism arises likely from single atoms or 
very small clustep-that are formed in the first stages of the preparation process and which do not give rise to further 
aggregation. Eiltja It was found that the fraction Xpara of Co which is deposited but is not detected by TEM increases 
as tco decreases. Accordingly, decreases as the amount of deposited Co decreases. On the other hand, the sample 
with tco = 1 iim has to,. = 0.71, that is, larger than 1/2, likely because of the predominant ferromagnetic coupling 
between particles. 

It is also remarkable that the average magnetic moment per atom for the whole sample, as obtained from Ms of 
Fig. 6, is smaller than the value for bulk Co {I.Tub per Co atom) for all samples and that it decreases as tco and, 
thus as. the average size of the particles decrease. This dependence is opposite ±q that observed for free Co clusters in 
beamaJ and also for Co particles of similar size supported in a solid matrixJja In those experiments, the measured 
magnetic moment per atom exceeded the bulk value and it was found to increase as the diameter of the particles 
decreases. 

The reduced value of Ms that we measure could be caused by an oxide layer at the surface of the particles, 
which orders antiferromagnetically. However, as we mentioned before, we did not find any trace of oxide in EELS 
measurements. Moreover, it is known that the exchange interaction between this layer and the magnetic core of the 
particles would also induce a net anisotropy on the latter .c3 This so-called exchange anisotropy leads to a shift of 
the hysteresis loops when the sample is cooled down in the presence of a magnetic field. For example, Peng and 
coworkerscEl have recently measured an exchange bias field as large as 10.2 kOe for CoO coated Co clusters having a 
diameter of 6 nm and 13 nm. By contrast, as we show in Fig. 4, the hysteresis loops measured after cooling the sample 
in zero field or in 5 T from room temperature are nearly identical, thus with no evidence for an antiferromagnetic 
order at the surface layer. Thus, we conclude that most of the particles are free from oxidation. 

It is however still possible that some of the Co atoms, in close contact with the AI2O3 matrix, have a weak chemical 
link with it. This chemical bonding can reduce the number of unpaired electrons and then quench the magnetic 
moment of the metal atom, as was shown by van Leeuwen and co-workers.B From our data, it is not possible to 
determine whether the atoms involved in the reduction of the average magnetic moment are located at the periphery 
of the particles or are those atoms which do not form aggregates, because the relative concentration of both increases 
as the average size of the clusters decreases. 

Therefore, in what follows, we approach the problem in a different way. As a starting point of the analysis we consider 
that the spheres have the bulk magnetization: Msb — 1.7/iB per Co atom, whereas the missing magnetic moment 
is exclusively attributed to the paramagnetic Co fraction. The contribution of the clusters to the net saturation 
magnetization of each sample then equals (1 — Xpara)Msb- It was obtained by subtracting the low-T paramagnetic 
tail from the total Mg- In this way, Xpara is also estimated. We list the results in Table I. The low-T paramagnetic 
magnetization is found to be compatible with a free spin 1/2 for all samples, which indicates that the isolated atoms 
have in average only one unpaired electron. 

In order to fit the magnetic data, it is necessary to know the fraction Xpara of paramagnetic atoms and their 
magnetic moment. 

B. Determination of the particle size distribution 

In this section we will try to determine the particle size distribution from the equilibrium magnetic properties of each 
sample and compare it with the results obtained by TEM. In Fig. 3 we have plotted the equilibrium magnetization M 
of two different samples having tco = 0.3 nm and 0.7 nm, respectively. We recall that, for a set of magnetic moments 
/i without anisotropy, M{H, T, fi) = MsL{ijH/T), where L denotes the Langevin function. If the anisotropy energy is 
taken into account, there is no analytical expression for M(H, T, fi) and the shape of the magnetization curve deviates 
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from tlie pure Langevin form when U /ksT is large. However, it can still be evaluated numerically, as was described 
m Ref.E3. As we said above, the good scaling on a single curve of data measured at different temperatures indicates 
therefore that the influence of the anisotropy is not very important. 

For a real sample, we have to average M(H,T, fi) over the appropriate distribution of particle's sizes. Comparing 
the calculated magnetization to the experimental data we shall try next to get information about this distribution. In 
order to directly compare our results with those obtained previously by TEM, we define g{D) as the distribution of 
number of particles having a diameter equal to D. For spherical particles /i — irAIsbD^ /6, where Mgb is, for the reasons 
given in the previous section, taken as equal to the saturation magnetization of bulk Co. We fit the experimental 
using the following expression: 

M{H, T) = Xpara^lB tauh j + (1 " XparajM^b- ^ g(D)V dD ^ ' 

taking for g a Gaussian distribution. 

For each sample, we fit only data measured at temperatures for which the two calculations, with and without 
anisotropy, give approximately the same result. We give an example of this method in the lower plot of Fig. 3 for a 
multilayer with tco — 0.3 nm. Above 30 K, the calculations performed with and without anisotropy almost coincide. 
For lower temperatures, close to Tg = 8.6 K, the experimental magnetization starts to deviate from the pure isotropic 
behavior, as happens for the data measured at T = 12 K. Even then, the experimental data are rather well reproduced 



by our calculations if we use the value of K^f / determined from the blocking of the ac susceptibility (see section IV C 
below). A list of the {D) values obtained from the fit for all samples is given in Table I. We find that {D) increases 
with tco, as expected. For most of the samples it is however larger by ten to forty percent than the values that were 
previously found by TEM. This discrepancy can be ascribed to the fact that the TEM experiments were performed on 
single layers deposited on a special carbon substrate, whereas we have measured multilayers prepared on a Si oxide. 
However, we have also measured a monolayer with tco — 0.7 nm and obtained almost the same magnetization results 
(see the upper plot of Fig. 3) as for a multilayer. An alternative explanation is that the saturation magneti: 
of the smallest particles is enhanced with respect to the bulk, as was found in similar systems of Co clusters 
However, even if we had used the maximum value of 2.3fj.B per Co atom, which was found by Respaud et al., (D) 
would have decreased by only a 10 percent, that is, within the uncertainty of the fitting procedure. In order to get 
the same diameters that were observed by TEM, Mg should be as large as 3/zb to 4/iB per Co atom for the smallest 
clusters. 

The fit of the magnetization curves is more sensitive to the value of {D) than to the width of the distribution a. In 
fact, it is possible to obtain a reasonably good fit by using a single Langevin curve with almost the same value of {D) 
(see Fig. 3, upper plot). In order to get a better estimation of cr, we have also fitted the equilibrium susceptibility, as 
obtained from the low field dc or ac measurements well above Tb, using the following expression 



/ giD) V {MlV/ikB T) dD 
JgiD)VdD 



^Ceg = r ^/n^TA^n ■ (3) 



This formula is valid also for pafticles with uniaxial magnetic anisotropy if, as it is the case for our samples, the 
anisotropy axes are not oriented.Ea It turns out that the equilibrium susceptibility is very sensitive to the presence 
of large particles in the distribution and therefore to a, as it is shown in the inset of Fig. 1. For this reason, the 
contribution of the paramagnetic moments to Xeq can be neglected in all cases. The values for {D) and a that are 
given in Table I are those which reproduce best both the equilibrium magnetization isotherms and the equilibrium 
susceptibility. The width of the distribution is found to be rather constant and in good agreement with the value 
found previously by, TEM. The slight increase of a as the average size of the aggregates decreases was also observed 
in the TEM dataJlJ In conclusion, the co-deposition of Co and AI2O3 gives us Co clusters of controllable size and 
with a narrow and nearly constant distribution of diameters. 



C. Magnetic anisotropy 



Next, we want to study the dependence of the effective anisotropy on the cluster size. The anisotropy of the particles 
can be estimated by comparing the average activation energy {U) to the average volume n/6{D)^. Usually {U) is 
estimated as (U) = 25kBTB, where Tb is the temperature of the maximum of the ZFC susceptibility. However, this 
pr^j^ure leads to an overestimation of the anisotropy because Tb depends not only on (U) but also increases with 
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In order to get more reliable values we need a way to obtain the full distribution of activation energies f{U), and 
then to find which value of U corresponds to particles with a diameter equal to (D). Fortunately, f{U) caB-|b€| 
directly determined from ac susceptibility data measured near the superparamagnetic blocking temperature Tb BO^Ell 
As mentioned above, the blocking occurs when the average relaxation time becomes of the order of l/w. It is therefore 
clear that the temperature dependence of x' and x" near Tb is determined by the distribution of U among the particles. 
In order to relate x' and x" to the distribution f{U) it is common to assume that those particles having U > Ub are 
fully blocked and the ones that fulfill the opposite condition arc in equilibrium. This hypothesis is reasonable because 
the relaxation time depends exponentially on U according to Eq. |l| For non-interacting particles, the relation that 
we were looking for reads as follows, 

/ XeqiU,T)f{U)dU+- x±{U,T)f{U)dU (4) 

Jo -J JUb 



X" ^^kBTxe,{T:Ut)f{Ut), (5) 

where Ub = fcBnn(l/cjro) is the activation energy of those particles having exactly t — t,, at a given temperature. 
Xeq — M^f^V/3kBT and X-ir-f -^sb/^^e// are the equihbrium susceptibility and the reversible (high frequency hmit) 
susceptibility, respectively.cS It follows from Eq. ^ that Uf{U) can be directly determined by plotting x" versus the 
scaling variable Ub- In Fig. 7 we show the result for a multilayer with tco = 0.7 nm. Similar results were obtained 
for the other samples. 

It is important to note here that f{U) is the fraction of the total magnetic volume occupied by particles having 
the activation energy U, since the susceptibility is mainly dominated by the contribution of the largest particles. 
Contrary, g{D) gives instead the number of particles of a given size. The two distributions are related as follows 

HU) = (6) 

' {dUldD)S^Vg{D)dD ^' 

For spherical particles U{D) ~ K^f f {tt/6) . Therefore f{U) is, apart from normalization factors, proportional 
to {U/Kef f^^^g [{GU/TrKeffY^^]- Using this relationship and taking as above a gaussian g{D), it is possible to fit 
x"{Ub)- The anisotropy constant is then simply the ratio between U{{D)) and V{{D)). Although the fit is rather 
good, we find that the function g that is extracted in this way from f{U) (or x") is systematically narrower than the 
size distribution obtained previously using the equilibrium magnetization and magnetic susceptibility. We will discuss 
later on the possible physical origin of this discrepancy. 

Before we comment on the variation of the anisotropy with the size, we would like to show that the distribution 
f{U) can also be obtained by a different method, which makes use of the ZFC and FC dc susceptibility curves 
measured at low enough magnetic fields. The difference between the ZFC and FC magnetization curves stems from 
the different contribution that the blocked particles make to each of them. Neglecting the weak variation of Adsb with 
T, this contribution only depends on T via the critical energy Ub which determines the relative number of blocked 
and superparamagnetic particles at a given temperature. Using the same approximation which led to Eqs. ^ and ^ 
for the ac susceptibility, it is possible to show that 

where Mi„ = — Mrev and Mrev is the magnetization brought by the reversible rotation of the magnetic moments. 
This expression is valid provided that the applied magnetic field is much smaller than the anisotropy field Hk = 
2Keff/Msb, as it is actually the case in our experiments. If this condition was not fulfilled, the activation energy 
would be a function of the field and of its orientation with respect to the easy axes of the particles. It is also possible 
to approximate Mirr{Ub) — Xeq{Ub)H. Therefore, Eq. gives an independent method to determine f{U). We plot 
in Fig. 6 the results obtained for an applied field of 10 Oe, which are in good agreement with the ac susceptibility 
data. In the same figure, data obtained for H = 100 Oe are also shown. In this case the maximum of the distribution 
shifts towards lower values of Ub, indicating that the activation energy decreases in a magnetic field. In addition, the 
distribution function broadens a bit as a result of the random orientation of the easy axes. 

We now come back to our main goal. The anisotropy constant is plotted in Fig. 8 as a function of the average 
diameter of the aggregates. It is interesting to compare these experimental data with the constant -K^e// that is 
estimated using only the intrinsic magnetocrystalline anisotropy of bulk Co. For hep Co, the stable phase for large 
particles, Keff equals the intrinsic uniaxial anisotropy constant K = 4.3 x 10^ erg/cm'^. However, the structural 
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characterization of all samples studied here shows that they crystallize in the fee phase. Therefore, we would expect 
that the intrinsic anisotropy of the particles in our samples would be smaller than for hep Co. For cubic anisotropyEEIEHI 
Keff = K/A, where K is the second order intrinsic anisotropy constant. Taking K — 2.8 x 10^ erg/cm^ for fee Co, 
this gives Keff — 7 x 10^ erg/cm^. Therefore, the values that we find for all samples are almost one to two orders 
of magnitude larger than expected for magnetocrystalline anisotropy. Furthermore, Keff is observed to increase as 
(D) decreases. The size dependence of the effective anisotropy follows approximately the following phenomenological 
expression 

Keff^K^ + ^ (8) 

with Koo — 5(2) X 10^ erg/cm^ and Kg = 3.3(5) x 10~^ erg/cm^. This result is robust in the sense that it does 
not change qualitatively if we use the average diameter found by TEM, instead of the values obtained from the 
magnetization data. The first term is close to K/A and can therefore be identified as the contribution of the intrinsic 
anisotropy. The second one is proportional to the fraction of atoms located at the periphery of the particles, which 
can be more than 80 % of all Co atoms for the smallest clusters studied here. Our experimental results indicate then 
that there exists a rather large contribution of the surface of the particles to the net anisotropy. 

The enhancement of the magnetic anisotttrtgi, of nanometer sized metallic particles with respect to the bulk has 
been previously reported by several authortfi'BEJ. For Co particles with diameters varying between 4.4 and 1.8 nm, 
Chen and coworkersQ obtained Keff which increases from about 5 x 10^ erg/cm^ to about 3 x 10^ erg/cm^. These 
values are even larger than ours. However, they are of the same order as the values that would have been obtained 
if we had used the tenmerature of the maximum of the ZFC susceptibility, as it was done by the authors. More 
recently, Respaud et al.Eil studied the anisotropy of Co particles of 1.5 and 1.9 nm by fitting the whole ZFC and FC 
magnetization curves, a method that can be considered as equivalent to ours. They found Keff ~ 8.3 x 10^ erg/cm^ 
and Keff — 7.3 x 10^ erg/cm'^, respectively, in reasonably good agreement with our data. The existence of a large 
surface anisotropy in metallic particles is thus well established experimentally. 

The origin of this extra anisotropjj-|lias-bcen related to the modification of the electrostatic and exchange interactions 
of the atoms located at the surface which depends largely on whether the surface is oxidized or not. However, 
as we argued above, the characterization of our samples by EXAFS, EELS, and XPS does not indicate the presence 
of an oxide layer .lla The same conclusion is derived from the hysteresis loops measured below Tb- Therefore, we 
have to consider how the properties of a "bare" metallic surface are modified with respect to the bulk. The valua 
of Kg that we have found is actually comparable to the perpendicular anisotropy measured in free Co surfaces 
It is commonly accepted-that this perpendicular anisotropy is related to the appearance of a large orbital magnetic 
moment on these atoms.LJ The 3d electrons become more localized at the surface and the localization gives rise to an 
increase of the orbital pHioment. The same theoretical interpretation can be applied to the atoms at the periphery of 
small metallic clustersD In this case, the enhanced anisotropy at the surface extends to the inner atoms via tho-Sitrong 
exchange interaction with them, which leads to an increase of the average anisotropy of even spherical clusters ElThis 
interpretation has been confirmed by X-ray magnetic dichroism experirpents performed on Au/Co/Au layersEj and 
more recently also on Co disk-like aggregates supported on Au surfaces.ll3. It was found that the orbital component 
171 L oi the total magnetic moment scales with the fraction of atoms located at the surface of the aggregates. For 
spherical clusters, as the ones studied here, we expect then that oc 1/{D), dependence that we have indeed 
observed for Kef / . We therefore conclude that the observed increase of Kef / is likely due to the increasingly localized 
character of the 3d electrons of the atoms located at the surface. 

Once the particle size distribution and the anisotropy are known, it is possible to predict the time-dependent 
magnetic response of the samples and compare it to the experiment. Examples of these calculations are compared to 
the experiments in Figs. 1, 3 and 5. The calculations account very well for the experimental data measured above Tb, 
as expected. They also reproduce in Fig. 1 the deviation of the FC susceptibility from the equilibrium susceptibility 
that takes place below 5 K. However, they reproduce neither the position nor the shape of the maximum of the ZFC 
susceptibility. Another example of this discrepancy is shown in Fig. 9, where we plot the experimental x' for a 
multilayer with tco — 0.3 nm and the values calculated (dotted line) with Eq. ^. Again, the width of the blocking 
transition is clearly overestimated by the calculations. 

We recall here that we have found that the activation barriers distribution is systematically narrower than the 
size distribution for all samples. As an example, in Fig. 10 we plot the size distribution g{D) of a multilayer with 
tco = 0.7 nm extracted from x" and directly observed by TEM. The horizontal scale for the former distribution is 
(&Ub/T^Kef fY^^ , with Keff = 10^ erg/cm"^. It is tempting now to attribute the "narrowing" of the blocking transition 
to the effect of the surface anisotropy. When K^jD 3> i^oo then U w KgS, where S = ttD^ is the surface of the 
particle. It follows then from Eq. || that f{U) oc D^g{D) and the width of the distribution of activation energies must 
then be smaller than when U oc V. Figure 10 shows indeed that when the same susceptibility data are represented 
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versus the variable {Ub/nKsy^"^ the ensuing size distribution is in better agreement with what it is found by TEM or 
from the equihbrium magnetization and susceptibiUty. In this way, we also obtain Kg which turns out to be between 
0.2 and 0.3 erg/cm^ for all samples. This value can be then used to recalculate the ac susceptibility and the ZFC 
magnetization. We find that the calculations performed with the same parameters a and (D) as before (see Table I) 
but taking U cx arc in much better agreement with the experiment (see Figs. 1 and 9). Although the width of 
the of size distribution is not always accurately determined, it seems that the influence of the surface anisotropy also 
modifies the shape of the susceptibility peak at the blocking. We conclude that the dynamical response of very small 
particles is therefore determined by the special physical properties of the atoms which are located at their surface. 

D. Influence of the number of layers: dipole-dipole interaction between the particles 

There has been some debate during the last years about the effect that r-the dipole-dipole interaction between 
magnetic nano-particles has on their relaxation times. Shtrikman and WolfarthcS and later Dormann et alxB predicted 
that the effective activation energy increases by an amount that depends on the number and spatial arrangement of 
the neighbor particles. By contrast, in the model proposed by S. M0rup and E. Tronca the interaction between the 
particles leads to a lower U. The experimental validation of one of these two models is complicated because, for some 
preparation methods, it is difficult to vary the density of particles in the sample without modifying the distribution 
of particle's sizes.Lj A different approach is to dissolve the particles in a fluid and to change the concentration by 
varying the amount of solvent. However, it is possible that the particles agglomerate in the-.fluid because of their 
mutual interaction, so that the interaction with the nearest neighbors is not greatly affected.ciJ 

The preparation method of our samples presents a number of advantages. We have seen that the average size can 
be controlled by changing the deposition time, but also the packing of the particles can be controlled. The TEM 
images show that the clusters in a layer do not agglomerate and, furthermore that tiie deposition of several layers of 
Co and AI2O3 leads to a self-organized spatial arrangement of the particles (see refJl3). For a multilayer each cluster 
has, in average, six nearest neighbors in the same plane, three above and another three below it. For tco = 0.7 nm, 
the average distance between nearest Co chisters in the same layer is A|| ~ 5.4 nm, whereas the distance to nearest 
neighbors in adjacent layers is A = 4.5 nmJij In this section, we compare the relaxation rate of two samples having 
both tco = 0.7 nm {{D) ~ 3 nm), but very different number of layers, namely 30 and only one. By going from a 
monolayer to a multilayer we certainly expect that the average energy of interaction of a particle with the others 
changes. The interaction energy between particles in adjacent layers is the largest and of the order of /x^/A'^ ~ 40 K. 
By contrast, in a sample with a single layer, each particles has, in average, only six neighbors coupled by a weaker 
interaction (/i^/A| « 20 K). 

In order to attribute any difference between the two samples to the effect of the interparticle interactions, it is very 
important to check beforehand that the sizes of the aggregates are the same in both. We showed in Fig. 3 that the 
equilibrium magnetization curves of the two samples are almost identical, and we compare in Fig. 11 the inverse 
of their ac susceptibility curves. Above Tb, the susceptibility follows the Curie- Weiss law, with identical values of 
C, which confirms that (D) and a are practically the same. By contrast, the Curie- Weiss temperature 9 is about 2 
times smaller for the monolayer, indicating that the average inter-particle interaction is notably reduced. It is also 
apparent that the blocking temperature of the monolayer is smaller than that of the multilayer. As we have done 
before, the activation energy of the two samples can be compared by plotting x" nieasured at different frequencies 
as a function of the scaling variable Ub- This comparison is shown in Fig. 12. The maximum of the curve for the 
monolayer is clearly shifted towards lower values of Ub with respect to the maximum obtained for the multilayer. Our 
data give strong evidence that the interaction between the aggregate layers tends to increase the activation energy 
of each particle, by an amount of about 200 K. This difference is of the same order of magnitude as the interaction 
energy with the six nearest neighbors in the multilayer. We also find that the relative width of [//([/) has the same 
value for the two samples, which confirms again that the distribution of particle's sizes is the same. 

E. Magnetic relaxation at low temperatures 

In the previous sections, the reversal of the magnetic moments has been treated as a classical process assisted 
by the interaction with a thermal bath. However, taken as a quantum variable, the spin of a magnetic particle 
S — MsV/gnB can in principle flip also by quant urp_tunneling across the barrier if the effective Hamiltonian contains 
terms which deviate from the uniaxial symmetry.Ej This possibility is very attractive because it would show the 
existence of quantum effects at the intermediate scale between the microscopic and the macroscopic worlds. Quantum 
relaxation can dominate over the thermal activation at very low temperatures, when the thermal population of the 
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first excited state doublet ±(5* — 1) becomes negligible, and should lead to a saturation of the relaxation rate to a 
nearly temperature i nde]3 | eiidq nt value. lj Such a saturation has indeed been observed in some systems of single domain 
particles in the past.ll3'lljl23'c3 

In this section, we present measurements of the relaxation of the remanent magnetization of an initially saturated 
sample. We have chosen the sample with the smallest Co clusters for two reasons; first, because the rate for quantum 
relaxation must be the largest for these clusters of only about 50 — 100 atoms; and second, because this sample shows 
the strongest anisotropy. The separation of the two lowest lying state doublets, which is roughly given by Qq ~ gfJ-sHk 
is then about 3 K, thus larger than the lowest temperature that our magnetometer can reach {Tmin — 1-7 K). We 
have measured the decay of the magnetization of the sample that takes place after a magnetic field of 5 T is switched 
off at different temperatures. The decay of Mr is approximately logarithmic in time. An important advantage of 
recording the relaxation at zero field is that it can then be easily calculated using our knowledge of the activation 
energies distribution. At zero field, the equilibrium magnetization is zero for all particles. Therefore, using the same 
approximation as before, the time dependent magnetization is given by 

M 1"°° 

M{t,T)^^ f{U)d{U) (9) 

where we have made the reasonable approximation that the magnetic moments of the particles are initially saturated 
by the magnetic field. The factor 1/2 arises from the reversible rotation of the magnetic moments for a random 
orientation of the easy axes, as in the Stoner-Wolhfarth modelrj As pointed out by Labarta et al.,E3 if the magnetic 
moments flip by a thermally activated process the relaxation curves measured at different temperatures should scale 
when plotted as a function of Ub- This plot also gives a picture of the relaxation at very long times, which are not 
experimentally accessible. Our experimental data, which we plot in Fig. 13 do indeed show a rather good scaling for 
the same tq = 10~^^ s that was obtained from the shift of the maximum of x" with frequency. The full line in the 
same figu re wa. s calculated with Eq. ^ using the distribution /([/) that we determined with the method described 
in section IV C . The scaling of the data confirms that the relaxation mechanism is classical (not tunneling) down to 
T = 1.7 K. 

In the inset of the same Fig. 13 we show the temperature dependence of the so-called magnetic viscosity 5*^, 
determined as the slope of the Mr vs ln(i) curves. Below about 2.5 K, Sr does not vary much with T. We note 
however that, according to Eq. ^, the magnetic viscosity is just 

and it is therefore proportional to f{U). The apparent saturation of S measured between 1.7 K and 2.5 K just reflects 
the shape of the distribution f{U), and it is indeed rather well described by the "classical" calculation. These data 
give aa example of how important it is to have information about f{U) in order to adequately interpret the relaxation 



data.i 



V. FINAL REMARKS AND CONCLUSIONS 



We have presented a detailed and extensive study of the magnetic properties of Co aggregates prepared by sequential 
deposition of Co and AI2O3. This preparation method enables us to control both the average size and the number 
of layers independently. We have shown that the distribution of activation energies can be accurately determined 
from ac susceptibility and ZFC-FC magnetization measurements. We have investigated the variation of the effective 
anisotropy as the size of the aggregates decreases from about 5 nm to below 1 nm. We find that K^ff scales with the 
fraction of atoms located at the periphery of the aggregates. The strength of the surface anisotropy is of the same order 
of what is found for free Co surfaces and we therefore attribute it to the increase of the orbital magnetic moment of 
these atoms. Furthermore, the activation energies distribution resembles the distribution of particle's surfaces rather 
than the volume distribution. For these small clusters, it is therefore more appropriate to write U — KgS than the 
"traditional" U = KeffV. Using the distributions of sizes and of activation energies that we have determined, we are 
able to give a quantitative account of all the equilibrium and time-dependent experimental quantities. We have also 
shown that the activation energy increases when the average number of nearest neighbors per particle increases, in 
agreement with the model of Dormann et al. Finally, the decay of the remanent magnetization of clusters containing 
only about 50 to 100 atoms is shown to proceed via a thermally activated mechanism down to the lowest temperatures 
investigated. 
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(a) (b) 

tco (nm) {D}" (nm) a" {D}^ (nm) a'' Xpara 

0.1 0.8(1) 0.35(5) 0.7 

0.2 0.83(20) 0.3 1.3(1) 0.3(1) 0.7 

0.3 1.4(3) 0.3 1.4(1) 0.32(5) 0.22 

0.4 1.4(3) 0.22 2.2(1) 0.2(1) 0.24 

0.7 2.9(6) 0.23 3.1(3) 0.2(1) 0.25 

1 4.2(8) 0.27 5.2(3) 0.25(5) 0.13 



TABLE I. Parameters of the gaussian distribution of particle's sizes obtained by TEM (a) and from the fit of the magneti- 
zation data (b). The width a of the distribution is given in units of the average diameter. The last column gives the estimated 
fraction of Co atoms which do not aggregate in particles. 
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Figure captions 

Figure 1. dc susceptibility of a multilayer with tco = 0.1 nm and N = 100 measured with a field of 100 Oe; open 
symbols, FC; closed symbols, ZFC. The lines represent the results of calculations performed with the parameters of the 
size distribution given in Table |: dashed line, equilibrium susceptibility; dotted line, ZFC susceptibility calculated 
taking U — K^ffV and K^ff = 2.4 x 10^ erg/cm"^; full lines, ZFC and FC susceptibilities calculated for surface 
anisotropy with Ks — 0.3 erg/cm^. Inset: Inverse suceptibility of three multilayers: (a), tco = 0.1 nm, N = 100; (b), 
tco = 0.3 nm, N = 40; (c), tco = 0.7 nm, N — 30. The lines represent the equilibrium susceptibility of (b) calculated 
for (D) = 1.4 nm and three values of a: 0.25 (upper curve), 0.3 (medium curve), and 0.35 (lower curve). 

Figure 2. Real and imaginary parts of the ac susceptibility of a multilayer with tco = 0.3 nm and N — 40. 

Figure 3. Equilibrium magnetization of multilayers with tco = 0.7 nm (a) and tco = 0.3 nm (b), measured at 
different temperatures. The lines represent the calculated results, (a): dotted line, pure Langevin curve for D — 3.1 
nm; full lines, results calculated averaging a Langevin curve over a Gaussian distribution of sizes (see Eq. ^) with 
a = 0.2 and three different values of the average diameter, (b): full line, as in the upper picture for a = 0.32 and 
(D) — 1.4 nm; dotted lines, equilibrium magnetization calculated for T = 12 K and T = 30 K with the same size 
distribution but for uniaxial anisotropy with U = irKgD^ and — 0.2 erg/cm^. 

Figure 4. Hysteresis loop of a multilayer with tco = 0.1 nm measured at T = 2 K after cooling the sample in zero 
field or in 5 T from room temperature. 

Figure 5. Left axis: Temperature dependence of the remanent magnetization and of the saturation magnetization 
(measured with H — 50 kOe) of a multilayer. The lines are calculated with Eqs. || and respectively using the 
parameters given in Table I and the distribution f{U) estimated from the blocking of the ac susceptibility. Right 
axis: Temperature dependence of the coercive field of the same sample. 

Figure 6. Variation of the reduced remanent magnetization (open symbols, right axis) and of the low tem- 
perature saturation magnetization Mg (closed symbols, left axis) with the amount of deposited Co for all samples 
studied. 

Figure 7. Imaginary part of the susceptibility of a multilayer with tco = 0.7 nm and = 30 plotted as a function 
of the scaling variable Us/kB — T In {1/ ujtq) , with tq = 10~^^ seconds. The full line is a fit according to Eq. ^ 
taking a Gaussian for g{D). Results obtained as explained in the text (cf Eq. ^) from ZFC-FC magnetization curves 
measured with two different magnetic fields are also shown for comparison. 

Figure 8. Size-dependence of the effective anisotropy constant for all samples investigated. The full line is a best 
squares fit of the data to Eq. |[ 

Figure 9. Real part of the susceptibility of a Co multilayer with tco = 0.3 nm and = 40 measured for two 
different frequencies. The dotted line is calculated for u)/2tt = 0.1 Hz with Eq. || using the parameters of Table | and 
taking U = irKeffD^/G, with Keff = 1.15 x lO'^ erg/cm^. The fuU lines are calculated taking U = irKgD^, with 
Ks = 2x 10-1 erg/cm^. 

Figure 10. The size distribution determined by TEM is compared to the distributions obtained from x" for 
two limiting cases where the anisotropy is either dominated by the intrinsic (volume) contribution (full line) or by 
the surface anisotropy (dotted line). The scaling in the horizontal axis gives respectively K^ff = 10^ erg/cm'^ and 
Ks = 0.33 erg/cm2. 

Figure 11. Inverse ac susceptibility of two samples with the same tco = 0.7 nm but different number of layers: 
open symbols, N — 30; full symbols, = 1. 

Figure 12. Scaling plot of x" for two samples with tco — 0.7 nm but different number of layers. For both 
To = 10"-^^ seconds. 

Figure 13. Time-dependent remanent magnetization of a (Co 0.1 nm AI2O3 3 nm)ioo multilayer plotted as a 
function of the scaling variable J/s/fcs = Tlnite/ro) with tq = 10^^^ seconds. The inset shows the temperature 
dependence of the magnetic viscosity. The full lines are calculated according to Eqs. and O. 
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